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Elvira R. SayfutyaroveQ 
Department of Chemistry, M.V. Lomonosov Moscow State University, Moscow 119991, Russia 

Alexei A. Buchachenkc0 

Department of Chemistry, M.V. Lomonosov Moscow State University, Moscow 119991, Russia and 
Institute of Problems of Chemical Physics RAS, Chernogolovka, Moscow District 142432, Russia 

Svetlana A. Yakovleva and Andrey K. Belyae\j^| 

Department of Theoretical Physics and Center for Advanced Study, 
Herzen University, St. Petersburg 191186, Russia 
(Dated: May 21, 2013) 

Charge-transfer cold Yb + + Rb collision dynamics is investigated theoretically using high-level ab 
initio potential energy curves, dipole moment functions and nonadiabatic coupling matrix elements. 
Within the scalar-relativistic approximation, the radiative transitions from the entrance ^4 1 E + to 
the ground X 1 E + state are found to be the only efficient charge-transfer pathway. The spin-orbit 
coupling does not open other efficient pathways, but alters the potential energy curves and the 
transition dipole moment for the A — X pair of states. The radiative, as well as the nonradiative, 
charge-transfer cross sections calculated within the 10 -3 — 10 cm -1 collision energy range exhibit 
all features of the Langevin ion-atom collision regime, including a rich structure associated with 
centrifugal barrier tunneling (orbiting) resonances. Theoretical rate coefficients for two Yb isotopes 
agree well with those measured by immersing Yb + ion in an ultracold Rb ensemble in a hybrid 
trap. Possible origins of discrepancy in the product distributions and relations to previously studied 
similar processes are discussed. 

PACS numbers: 34.20.Cf,34.70.+e,34.50.Cx 



I. INTRODUCTION 

Ongoing efforts in the creation of cold atomic and ionic 
ensembles, as well as in their subtle manipulation, have 
been recently merged in preparation and study of the 
hybrid ion-atom systems [TMlO). This signifies the im- 
portant step towards probing, understanding and explor- 
ing ion-neutral interactions and collisions at sub-Kelvin 
temperatures. Elastic and momentum transfer collisions 
may serve as the means of cooling of both neutral [TTITB] 
and charged [3HH[S] components. Ion-atom collisions re- 
sulting in charge transfer are of interest for a prototype 
of reactive processes and suggested to lead to gas-phase 
"metal-to-insulator" transitions at microKelvin temper- 
atures [15] . 

Elastic, as well as resonant charge transfer, ion-atom 
collision dynamics is well understood theoretically [13, 
116] . High-energy limit mediated by a short-range ex- 
change interaction transforms to the so-called Langevin 
regime when a collision energy becomes comparable to 
a long-range induction interaction. When a collision 
energy decreases further, the quantum regime with its 
threshold laws for each partial wave is attained. This 
limit has been thoroughly considered within the mul- 
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tichannel quantum defect theory fPTH2"0"] , but not yet 
reached experimentally. 

Preparation of cold ion-atom systems is achieved in a 
hybrid trap by bringing in contact atomic and ionic en- 
sembles or by creating ions via photoionization of cold 
atoms [11113 0113 [II [EE]- Collision- or photo-induced 
processes in Coulomb crystals can also be studied below 
10 K [I] [5H H3] providing that the neutrals are cooled 
or decelerated [Mj. However, in all cases the motion of 
ion(s) in the trapping field creates an unavoidable inher- 
ent source of a kinetic energy that imposes the limits on 
cooling [25 . Effective temperatures of a few tens of mil- 
liKelvin attained experimentally [U [5] still correspond to 
the Langevin regime with its complex dynamics sensitive 
to the global ion-atom potentials, transition moments 
and coupling matrix elements. In addition, laser cool- 
ing creates appreciable amount of electronically excited 
species, which further complicate the dynamics [7[ 126]. 

This explains the need in thorough theoretical studies 
that account for specific features of the particular per- 
spective system. Resonant charge transfer has received 
more attention, with the focus on the model alkali or 
alkaline-earth (e.g., Refs. QH H3 [5S]) and Yb+ + Yb 
[29) systems. For nonresonant case, few combinations of 
alkaline-earth ions and alkali atoms [13j [HI f30T433J . as 
well as Yb + , Ba + + Ca collisions [5J [26], were consid- 
ered. Accurate ab initio calculations that back most of 
these studies revealed the complexity of charge-transfer 
pathways in the above systems. Strongly bound excited 
electronic states may well be involved in the dynamics 
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by the nonadiabatic and/or spin-orbit couplings, or even 
directly upon laser cooling, preventing the use of simple 
two-state models [5] [S5J EHl IMJ 133] • Collision dynamics is 
complex itself being dominated by the resonances typical 
to the multiple-partial wave Langcvin regime. 

In the present paper we expand the theoretical expe- 
rience considering the charge transfer in cold collisions 
between the ground-state Yb + ion and Rb atom. Ex- 
perimental studies by Kohl and coworkers [3J [3] , further 
extended to excited Yb + ions [8] , demonstrated the pos- 
sibility of cooling an Yb + ion in the ultracold Rb en- 
vironment, gave the charge-transfer rate coefficient of 
the order of 10 -14 cm 3 /s (the smallest one measured so 
far below 1 K) and characterized the product distribu- 
tions. This system is attractive for further experimental 
studies since both Yb( 1 S) and Rb( 2 S) neutrals can be 
brought to degenerate gases [34TI35] and combined with 
Rb + ( 1 S) and Yb + ( 2 S) ions for studying charge trans- 
fer from both "sides". Rich isotope variety enables one 
to address the isotope and hyperfine structure effects. 
We performed accurate ab initio calculations of the po- 
tential energy curves, dipole moments and nonadiabatic 
coupling matrix elements for the lowest electronic states 
of the (YbRb) + ion. We showed that within the frame- 
work of the scalar-relativistic approximation the charge 
transfer can be considered by means of the two-channel 
model and that its radiative pathway is barely dominant. 
The vectorial spin-orbit interaction affects the electronic 
structure of the (YbRb) + ion, but does not produce a re- 
markable effect on the low-energy radiative charge trans- 
fer. Quantum scattering calculations within the range of 
collision energies relevant to the experiment show a rich 
resonance structure associated with orbiting resonances 
typical for the Langevin regime. Calculated rates agree 
well with the measured ones [3], though reveal less pro- 
nounced isotope effect. On the other hand, our study 
of single-collision dynamics results in some deviation in 
interpretation of the product distributions deduced ex- 
perimentally. 



II. ELECTRONIC STRUCTURE 

The electronic structure of the (YbRb) + ion was in- 
vestigated ab initio using the MOLPRO program pack- 
age [37]. All calculations were performed within the C2 V 
symmetry group with the origin of the electronic coordi- 
nates placed at the center of nuclear mass of the system 
composed of 174 Yb and 85 Rb isotopes and z axis oriented 
along the internuclear vector R. 



A. Techniques 

Scalar-relativistic (SR) calculations were performed us- 
ing the small-core 28-electron relativistic core potentials 
ECP28MDF with the supplementary contracted basis 
sets for both Yb [35] and Rb [35]. The sets of dif- 



fuse primitives of each of the spdf types with the ex- 
ponents continuing two lowest exponents of the stan- 
dard bases as an even-tempered sequence were added 
at each center to improve the description of induction 
and dispersion interactions. To evaluate the potential 
energy curves (PECs) of the states that are the lowest 
ones in their symmetry (spatial and spin) representation, 
the restricted version of the coupled cluster method with 
singles, doubles and noniterative triples, CCSD(T), was 
employed with the restricted Hartree-Fock (RHF) ref- 
erence. This type of the calculations always included 
the set of 3s3p2d2flg bond functions (bf) [30] placed 
at the middle of the internuclear distance R and coun- 
terpoise correction for the basis set superposition error 
[IT] . The Yb(4s 2 4p 6 4d 10 ) shells were included in the 
core, whereas all the rest electrons were correlated ex- 
plicitly. The lowest excited singlet states were calculated 
by means of the equation-of-motion approach in coupled 
clusters with singles and doubles (EOM-CCSD) imple- 
mented in the same way. For excited triplets the mul- 
tireference configuration interaction (MRCI) calculations 
with the Davidson correction |42j were performed with 
the reference wave functions built by the state-averaged 
complete active space multiconfigurational self-consistent 
field (CASSCF) method with active space spanned by the 
Yb(6s6p) and Rb(5s) atomic orbitals. The core set in the 
MRCI calculations consisted of Yb(4s 2 4p 6 4d 10 5s 2 5p 6 ) 
and Rb(4s 2 ) atomic orbitals, whereas the 4/ 14 shell of 
an Yb atom was correlated as fully occupied. 

The similar MRCI method was used for spin-orbit (SO) 
calculations, but with the ECP28MWB effective core po- 
tential for Yb containing the SO part 43! and supplemen- 
tary segmented basis set [H] augmented by the s2pdfg 
diffuse functions [45] . 

The calculations were performed on the fine grid of the 
internuclear separation R from 1.9 to 40 A. 



B. Scalar-relativistic results 

The results of SR CCSD(T), EOM-CCSD and MRCI 
calculations are shown in Fig. 1(a). At low collision en- 
ergies it is sufficient to consider the states correspond- 
ing to the three lowest dissociation limits: (i) Yb( 1 S) + 
Rb + ( 1 S) that represents the final charge transfer (CT) 
channel, (ii) Yb + ( 2 S) + Rb( 2 S), the initial channel, and 
(iii) closed Yb*( 3 P°) + Rb+( 1 S) CT channel that lies 
slightly above the entrance. Our best estimations for en- 
ergies of excited limits 16279 cm" 1 (EOM-CCSD) and 
18646 cm -1 (MRCI) correspond well to the centers of 
the measured fine-structure multiplets, 16750 and 18869 
cm -1 , respectively [35]. It is evident that direct CT can 
only occur through - Y 1 E + interactions, either 

nonadiabatic or dipole ones. 

To obtain the most accurate SR PECs for the rel- 
evant X 1 ^, a 3 S+, 6 3 n and 2 3 £+ states, we 
took CCSD(T)/bf results for X, a and b as the refer- 
ence. The A-state PEC was obtained by adding EOM- 
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TABLE I: Parameters of the lowest PECs of the (YbRb) + 
ion: Equilibrium distance, R e , state binding energy, D e , lead- 
ing LR coefficient, C n , energy of the minimum with respect to 
the ground-state asymptotic limit, U e , and state dissociation 
asymptote, Una, with respect to the ground asymptotic limit. 



FIG. 1: (color online) Potential energy curves and dipole mo- 
ment functions for the lowest electronic states of the (YbRb) + 
ion. (a) Scalar-relativistic PECs. In addition to the states cor- 
relating to the three lowest asymptotic limits discussed in the 
text, 3 1 E + and 1 II states correlating to the fourth Yb( 1 P°) 
+ Rb+( 1 S) asymptote (26172 cm* 1 , cf. 25068 cm -1 [46]) are 
shown, (b) PECs for the lowest SO-coupled 0* and 1 states. 

(c) The scalar-relativistic A — X transition dipole moment. 

(d) The SO-coupled A - X, al - X and 3 0+ - X transition 
dipole moments. 



CCSD excitation energies to the ground PEC, whereas 
the MRCI results for II — E splitting were used to ob- 
tain 2 3 E+ PEC from the CCSD(T)/bf fe 3 II one. These 
ad hoc procedures aim to compensate the lower accu- 
racy of EOM-CCSD and MRCI methods and incorporate 
the basis set superposition error correction. The result- 
ing points were interpolated by cubic splines and joined 
smoothly to an analytical long-range (LR) function rep- 
resented by the lowest-order —C n /R n asymptotic term. 
Finally, the PECs were shifted in energy to reproduce the 
true asymptotic limits. The parameters of so-obtained 
PECs are presented in Table [I] 

Accuracy of the ab initio PECs can be indirectly ver- 
ified at long distances. LR behavior of the X and a, 
A potentials is dominated by Yb and Rb induction, re- 
spectively. The finite-field CCSD(T) calculations esti- 
mated the static dipole polarizabilities of neutral atoms 
as «Yb = 142.2 and amj = 318.2 a.u., in perfect agree- 
ment with the recommended values ayb = 139±7 [4TIFI8] 
and a Rb = 318.8 ± 1.4 a.u. USED]. The leading induc- 
tion coefficients C4 — ax/2 are equal to 71 and 159 a.u., 
respectively, in reasonable correspondence with the fits 
to ah initio PECs, see Table [T] Asymptotic dependence 
of the & 3 LT and 2 3 E + PECs originates from the charge- 
quadrupole Rb + + Yb*( 3 P°) interaction. The results 
of Ref. [48] obtained within the more accurate ab initio 
approach allowed us to estimate the corresponding C3 
coefficients as -7.2 a.u. and 14.4 a.u., respectively. They 
agree with the fitted coefficients only qualitatively, but 
the LR behavior of the states correlating to the third 
dissociation limit is not crucial in the present context. 
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b Right minimum, corresponds to the A 1 E + SR state. 
c Left minimum, corresponds to the 6 3 II SR state. 



d Avoided crossing maximum. 



The permanent and transition dipole moments for the 
pair of 1 E + states were calculated by EOM-CCSD and 
MRCI methods, while the finite-field CCSD (T)/bf calcu- 
lation was possible also for the ground X state. All the 
methods give similar results, though the MRCI A-state 
moment reveals the signatures of the mixing with higher 
lying states at the distances shorter than 4 A. The A — X 
electronic transition dipole moment oIax responsible for 
radiative CT is shown in Fig. 1(c). 

In addition to the transition dipole moment, two 1 E + 
states are coupled by the nonadiabatic coupling matrix 
elements (NACMEs) as defined in Sec. Ill The first-order 
radial NACME was computed using the two-point finite- 
difference procedure [37] using the MRCI vectors built 
on the CASSCF wave functions obtained with averag- 
ing over the two lowest 1 E + states. It was checked that 
neither inclusion of the third state of the same symme- 
try into CASSCF averaging nor the use of more accurate 
three-point differentiation procedure alter the results re- 
markably. The second-order radial NACME was approx- 
imated as the derivative of the first-order one [5T]. For 
two states of the same symmetry the first-order angu- 
lar NACME, proportional to the orbital electronic angu- 
lar momentum operators L x and L y , vanishes by parity. 
Non-vanishing is the second-order NACME represented 
by the l? x + L 2 y operator. The one-electron part of the l? y 
operator (L x is the same by symmetry) was calculated as 
the matrix element on the CASSCF wave functions. As 
shown in Fig. 2, the first-order radial coupling is expect- 
edly much stronger than the second-order and angular 
ones. The mixing with the higher lying electronic states 
at short distances already noted for dipole moments af- 
fects NACMEs as well. 
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the corresponding dipole moment is negligible. 

III. NONRADIATIVE CHARGE TRANSFER 

The nonradiative charge transfer process in low-energy 
Yb + + Rb collisions can occur due to nonadiabatic tran- 
sitions between the initial and the final molecular states 
induced by nuclear motion. We estimated its probabil- 
ity using the SR ab initio picture presented above, see 
Fig. 1(a). In this approximation, it is sufficient to con- 
sider only two states: the initial A and the final X ones, 
both of the 1 S + (A = 0) symmetry, whereas the spin 
restriction forbids the low-energy CT in the triplet man- 
ifold. 



FIG. 2: (color online) Nonadiabatic couplings between 
the X X E + and v4 1 E+ states: the first-order radial NACME 
(dashed line), the second-order radial one (dotted line) and 
the second-order angular NACME multiplied by 50 (solid 
line). 



C. Spin-orbit coupling 

The state-interacting MRCI SO calculations [52] were 
performed in the space spanned by the scalar-relativistic 
states correlating to three lowest dissociation limits. The 
SO matrix elements were evaluated for Breit-Pauli oper- 
ator using the inner part of the MRCI wave functions. 
The PECs described above were used as the diagonal SR 
part of the full Hamiltonian matrix. Relevant SO-coupled 
PECs are shown in Fig. 1(b) and characterized in Table 
III 

The states correlating to the lowest dissociation lim- 
its are weakly affected by the SO coupling. The only 
qualitative effect, the splitting of the a 3 £ + state into the 
a0~ and al components, is not significant. Much more 
pronounced is the interaction with the higher states cor- 
relating to an excited Yb* atom. The asymptotic limit 
splits into three components that correspond to the 3 P°, 
j = 0, 1 and 2, Yb* fine-structure levels. The lowest 
3 Pq one gives the single 0~ component that may interact 
with the a0~ state at short distances, whereas the sec- 
ond 3 P° gives + and 1 components. The former trans- 
forms 6 3 II — A 1 Y, + crossing into the avoided one giving 
the double-well A0 + PEC. The latter perturbs the repul- 
sive wall of the al potential. These changes also affect 
the transition dipole moments. The A — X transition 
moment is similar to its scalar-relativistic precursor at 
distances larger than the avoided crossing point R ~ 6 A 
and rapidly falls down at shorter distances, see Fig. 1(d). 
Instead, a significant dipole coupling occurs with the up- 
per 3 + state, which acquires charge-transfer character 
as the result of the crossing. The sum of A — X and 
3 + — X transition moments is therefore similar to the 
scalar-relativistic A 1 S + one. The SO interaction 

also allows the transition from the al to the X state, but 



A. Theory 

The nonadiabatic nuclear dynamics was studied 
within the formalism of the standard adiabatic (Born- 
Oppenheimer) approach described, for example, in 
Refs. [53Tl56| . The approach is based on a fundamen- 
tal simplification, the Born-Oppenheimer separation of 
electronic and nuclear motions, leading to a fixed-nuclei 
electronic structure calculation and an appropriate treat- 
ment of nuclear motion based on the data calculated in 
the first step. This separation results in the total wave 
function expanded in terms of products of the electronic 
fixed-nuclei wave functions $ j(r, R), the angular nuclear 
wave functions and the radial nuclear wave functions, r 
being a set of electronic coordinates. The electronic wave 
functions are calculated as described in the previous sec- 
tion. The angular nuclear wave functions are expressed 
via the generalized spherical harmonics [55, 55] . The ra- 
dial nuclear wave functions obey the system of coupled 
channel equations [55J [55] . In the present case of low- 
energy Yb + + Rb collisions, when nonadiabatic transi- 
tions occur only between two 1 E + states, the coupled 
channel equations for the radial nuclear wave functions 
F-' E (R) are reduced to the following equations (in a.u.) 



1 d 2 
2n dR 2 



Uj(R) 



J(J+1) 
2fiR 2 



E 



pj,E = 



dR' 



I** 



dF, 



J,E 



dR 



1 

2/z 



k 



'clR 2 ' 



Wk)F 



J.E 



J,E 
k 



Here E is the collision energy measured from the asymp- 
totic limit of the initial state J 4 1 E + , J is the total an- 
gular momentum quantum number, which in case of the 
1 S + states represents simultaneously orbital and rota- 
tional momenta of the nuclei, [i is the reduced mass of 
the nuclei, j and k indexes runs over A and X, Uj is 
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the adiabatic PEC and three terms in the right-hand 
side contains first- and second-order radial and angular 
NACMEs, respectively, defined as the integrals over the 
electronic coordinates r and calculated as described in 
the previous Section. 

The coupled channel equations have their simplest 
and standard form due to the choice of coordinates: the 
Jacobi coordinates in which the vector R connects the 
nuclei (for a fixed-nuclei treatment) and the set of elec- 
tron coordinates r is defined from the center of nuclear 
mass, the coordinates employed in the ab initio calcu- 
lations described above. In these coordinates a special 
care, e.g., by means of the reprojection method [5rjVl58j . 
should be taken in the asymptotic (R — > oo) region for 
calculating of nonadiabatic transition probabilities due 
to the fact that some radial NACMEs may have non- 
vanishing values, but in the present case all the treated 
NACMEs have zero asymptotes, see Fig. 2. 

Due to the large energy splitting and small values of 
NACMEs, the nonadiabatic transition probabilities are 
expected to be small at low collision energies, and further 
simplification can be achieved by using the perturbation 
theory. The probabilities for a nonadiabatic transition 
4^1 due to the radial NACMEs (of both orders) can 
be approximated as j59] 



J,E 



vJ,e dF x 



F J, E dF^ E 
tx dR 



dR 



(2) 



F^ A AR) being an unperturbed scattering radial wave 
function (distorted wave) in the channel X (A) normal- 
ized by the probability current. For transitions induced 
by the angular NACME of the second order, the equa- 
tions are similar to the coupled channel equations in a 
diabatic representation, and the perturbation theory pro- 
vides the following nonadiabatic transition probability 



B. Cross sections 

The calculated A — > X nonradiative cross sections are 
presented in Fig. 3. All calculations were performed for 
the 85 Rb isotope. The 174 Yb isotope is assumed unless in- 
dicated explicitly. It is seen that the CT process induced 
by the angular NACME is roughly 14 orders of magnitude 
less efficient than that induced by the radial couplings. 
The reason is that the total electronic orbital momenta 
L, though not the good quantum numbers, are zero in 
the asymptotic R — > oo limit. (In the united ion limit, 
Bh + , at least the X state should also have L = |60|.) 
Nonvanishing angular NACMEs can only arise from the 
minor admixtures of excited states with different angular 
structure to the adiabatic electronic wave functions. In 
turn, the cross section due to the radial NACMEs is small 
on its own: the corresponding rate coefficient is roughly 
of the order of 10 -26 cm 3 /s at 1 K temperature. This is 
the consequence of both the large adiabatic splitting and 
weak interaction of the A 1 Y 1 + and X 1 T, + states, resulting 
in small values of the radial NACMEs. 
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The nonradiative CT cross sections and tr^ G due 

to the radial and the angular NACME are computed as a 
sum over the total angular momentum quantum number 
J: 



FIG. 3: (color online) Nonradiative CT cross sections for 
Yb + + Rb collisions induced by NACMEs of different ori- 
gins. The uppermost (blue) trace shows the cross section 
for 172 Yb+ + Rb collisions due to the radial NACMEs. The 
second from the top (black) trace represents the same cross 
section for 174 Yb + + Rb collisions (multiplied by 10 -3 ). The 
bottom (red) trace depicts the cross section for 174 Yb+ + Rb 
collisions due to the angular NACME (multiplied by 10 8 ). 
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gi being a statistical probability for population of the ini- 
tial channel i. For the entrance A 1 Y, + channel, = 1/4. 
Because of low collision energies treated the radial wave 
functions have to be calculated up to large internuclear 
distances, up to several thousands Angstroms in the 
present case. 



The nonradiative cross sections exhibit numerous reso- 
nances whose positions in the radial and angular counter- 
parts are identical for the same Yb isotope. This struc- 
ture, typical to Langevin regime [mi ITt)I |2"5H3T)1 13"3"] . arises 
from the tunneling through centrifugal barriers of the ef- 
fective potential energy in combination with the polar- 
ization potential for the initial channel, as is shown in 
the next Section. Isotope substitution alters the cross 
section insignificantly. 
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To conclude, the large mass of the system, vanishing 
electronic angular momenta of the X and A states and 
their large energy separation in the interaction region 
suppress the nonadiabatic couplings and make the direct 
nonradiative CT process very inefficient. Inclusion of the 
SO couplings does not change the conclusion. For this 
reason, we did not exploited more thorough treatments 
of the nonadiabatic dynamics. 



IV. RADIATIVE CHARGE TRANSFER 

The dipole coupling between A and X states should 
provide more plausible radiative CT mechanism. As far 
as the nonadiabatic effects are negligible, we estimated its 
efficiency using the standard adiabatic picture and com- 
monly accepted quantum scattering approaches. Follow- 
ing the previous works on cold ion-atom collisions, we did 
not consider the effects of external fields and hyperfine 
structure. 



A. Theory 

Detailed description of the collision-induced sponta- 
neous radiative transition from electronic state % to elec- 
tronic state / is given by the frequency-resolved cross sec- 
tion da R (E) /duj, where the collision energy E and the fi- 
nal energy E' define the emission frequency hw = E — E' . 
Depending on the final PEC and the final energy, tran- 
sitions may occur into the continuum (charge transfer 
in the strict sense, hereafter "charge exchange", CE) 
or bound (radiative association, RA) levels. Within 
the Fermi Golden Rule (FGR) approximation the cor- 
responding cross sections are expressed as 



da CE {E) _ 8ir 2 



did 



3 fc 3 fc' 



J 



(2J+1)J2h, 



) 2 

'EJ,E'J' 



and 



da HA (E) 
duj 

-EE 



s^ 2 

3 



fc 3 



9iOt 



3 ^(2 J 

J 



,8{u - w E ,v>j')- 



(5) 



(6) 



Here J and J' are the nuclear angular momenta in 
the initial and final states, v' is the vibrational quan- 
tum number of the final state with the energy E v iji, 
Hwe,v>J' = E — E v rji, a is the fine structure constant, 
k 2 = 2(iE, k' 2 = 2[iE', S(x) is the Dirac function and 
Hjji is the standard Hdnl-London factor. 

The transition dipole moment matrix elements on the 
radial nuclear wave functions are defined as 



Dej,E'J> 
Dej,v'J> 



/>oo 

jf F?< E (R)d if (R)Ff> E '(R)dR, 

poo 

J F: 7 - E (R)d lf (R)Ff-'(R)dR. 



The choice of the coefficients in Eqs.|5]), (|6| implies the 
normalization of distorted waves , F/ ' a used m 
Ref. [53]. 

The total radiative cross section is given by 
a§(E) = af f E (E) + a^ f A (E) 



dco 



da?f(E) daf/(E) 



did 



"o daf f E (E) 

duo 

duj 



da? f A (E) 
dui 



(7) 



where w max is the maximum frequency that corresponds 
to a transition to the ground level v' = 0, J' = 0, whereas 
Wo frequency does to the ground-state dissociation limit. 
After substitution of Eq. ^ into Q , integration over u> 
with (5-function transforms the last term of the latter into 
a sum over the discrete J' quantum numbers. 

The closed expression for the total elastic cross section 
follows from the optical potential (OP) approach (see, 
e.g., Refs. [64 1I66H68] and references therein). Combined 
with the distorted wave approximation, it reads [64] 

= Ti9i E( 2J + ^i 1 - expf-^^)]} (8) 



with 



2?r h 3 

r 

3 k 



/ r u d:,:R) l , I i 4 i ' 1 <IH. (9) 
Jo 



Vj( e ) = -■ ■■■<'■ 

For the - or AQ+-XQ+ transitions under 

consideration, Hjji — J(2J+1) _1 if J' = J — 1, Hjj> = 
(J + 1)(2J + l)" 1 if J' = J + 1 and Hjj. = otherwise. 



B. Cross sections 

The radiative cross sections defined above were com- 
puted numerically using fine radial and energy grids and 
strict convergence criteria for partial wave summations. 
Accuracy of the OP cross sections was estimated as 
3%, whereas the accuracy of the FGR calculations were 
within 5-8% due to additional errors in integration over 
lo and convergence of the wave functions very close to the 
dissociation limit of the ground state. 

The total OP radiative cross section for the SR model 
(all the parameters, A 1 S + , A" 1 E + PECs and dAx tran- 
sition moment function, are taken from the SR calcu- 
lations) for the collision energy range 0.001 - 10 cm -1 
are shown in Fig. 4(a) as the upper trace (multiplied by 
two). It gradually declines with collision energy increase 
in agreement with the Langevin capture cross section 
ul = 2-KyJC^/E being five orders of magnitude smaller. 
The radiative CT cross section exceeds the nonradiative 
ones by 14 orders of magnitude and bears similar reso- 
nance structure. Analysis of the individual partial wave 
contributions allowed us to correlate the resonance ener- 
gies with the positions of centrifugal barriers at certain J, 
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as is shown in Fig. 5 for the low-energy cross section for 
collisions involving 172 Yb and 174 Yb isotopes. The cross- 
section structure consists of strong narrow resonances, 
which correspond to the tunneling deep under the cen- 
trifugal barriers, and weak broad resonances slightly be- 
low or above the barriers. All these findings signify the 
Langevin character of the Yb + + Rb cold inelastic colli- 
sions. 





■—-^■-•0,1 (L 
(b) 



1E-3 0.01 0.1 1 10 
E(cm"') 



36 38 40 42 44 
E(cm" 1 ) 



FIG. 4: (color online) Total radiative A — X cross sections 
for the Yb + + Rb collisions calculated within the OP ap- 
proximation, (a) Low-energy results for SR (upper trace, 
multiplied by two) and SO (lower trace) models. Asterisks 
represent the RA (free-bound) cross section from the FGR 
calculations for the SO model. Dashed line corresponds to 
the Langevin inelastc cross section multiplied by 10 -5 . (b) 
Cross sections near the potential barrier top of the A0 + state 
(40 cm -1 ). Dotted and solid lines - SR and SO models, re- 
spectively. Shaded area presents the results of the "mixed" 
model calculations with the SO-coupled PECs and SR transi- 
tion dipole moment. Asterisks indicate the resonance features 
of the SO cross section associated with the potential barrier. 



The total OP cross section computed within the SO 
model the SO-coupled X0+ , A0+ PECs and the d AX 
transition moment) are also presented in Fig. 4(a) (lower 
trace). Test calculations showed that the al — X0 + radia- 
tive transition has a negligible probability. The SO model 
gives the cross section by 10% larger than the SR one with 
the same energy dependence, except the "bumps" at ca. 
0.001 and 0.2 cm -1 . Resonance structures of two traces 
are not identical, but the correlations are clearly visible 
at least below 0.01 cm -1 and were confirmed selectively 
by the partial wave analysis. However, the latter was of 
little help for seeking the effect of the potential barrier 
separating two wells of the SO-coupled A-state PEC, see 
Fig. 1(b). 

To elucidate it, we performed OP calculations with 
the "mixed" model that uses SO-coupled PECs but SR 
transition dipole moment. The results of three models 
are compared in Fig. 4(b) for collision energies close to 



E 

o 



: 19 120121 



27123291 3D3132 

i 




E(cm" 1 ) 



0.005 0.010 



FIG. 5: (color online) Low-energy resonance structure of the 
total radiative A — X cross sections (the SR model, OP cal- 
culations). Vertical lines indicate the positions of centrifugal 
barrier top at J values given at the top of each plot. The res- 
onances are assigned to J by arrows or J values given at the 
bottom of each plot, (a) 172 Yb+ + Rb collisions; (b) 174 Yb+ 
+ Rb collisions. 



the potential barrier top (40 cm 1 , see Table Pre- 
serving the Langevin-type structure similar to that of 
the SR cross section, the "mixed" cross section is dom- 
inated by additional very intense features significantly 
broadened at the background. They were attributed to 
the potential barrier tunneling or overbarrier resonances, 
since the penetration into the short-range well enhances 
the radiative transition probability (both Franck-Condon 
overlaps and transition moment value are more favorable 
here than in the long-range well, see Fig. 1). It is not the 
case for the consistent SO model, in which the transition 
moment decline suppresses the radiative transitions from 
the short-range well. Some features associated to the 
potential barrier, however, survive in the SO cross sec- 
tions (most visible are marked by asterisks in the figure) . 
The same "mixed" model was used to confirm that the 
"bumps" mentioned above are also due to the potential 
barrier. As additional comments to Fig. 4(a), increase of 
the collision energy above the potential barrier does not 
lead to any drastic change of the cross section because 



the centrifugal term smoothes it for many partial waves 
and increase of the SO cross section with respect to SR 
one reflects the variation of the PECs rather than of the 
transition moment. 

The results of the FGR calculations performed for the 
SO model at selected collision energies are presented in 
Table [TTJ The total FGR cross section agrees with the 
OP one within the numerical accuracy justifying the op- 
tical potential approximation. Radiative association pre- 
vails over charge exchange: within the energy range un- 
der study, transitions to the bound rovibrational levels 
of the (YbRb) + ion amount ca. 70% of the total radia- 
tive charge transfer, see Fig. 4(a). Figure 6 compares 
A — X spontaneous emission spectra at the collision en- 
ergies 0.01 and 1 cm -1 . Increase of the collision energy 
leads mostly to "rotational" congestion of the spectrum 
due to involvement of higher partial waves. It is worthy 
of noting that all transitions are concentrated in the rela- 
tively narrow frequency range, from 500 cm -1 above the 
ground dissociation limit to 2000 cm -1 below (ca. 12% 
of the available frequency range from to 20250 cm -1 ). 



TABLE II: Radiative cross sections (in 10 -17 cm 2 ) for charge 
transfer in the Yb + + Rb collisions calculated within the SO 
model at several collision energies E. Percentage contribu- 
tions of RA and CE processes are given in parentheses. 



E, coT 1 




FGR 




OP 


RA 


CE 


Total 


Total 


0.01 


10.2 (71) 


4.2 (29) 


14.3 


15.2 


0.1 


3.0 (69) 


1.4 (31) 


4.3 


4.5 


1 


0.9 (69) 


0.4 (31) 


1.3 


1.3 


3 


0.05 (67) 


0.03 (33) 


0.08 


0.08 



V. DISCUSSION 

In this Section, we first compare our theoretical results 
with experimental data available for Yb + + Rb collisions. 
Then we discuss them in a broader context of a few previ- 
ous studies, mainly theoretical, on the cold CT processes 
in analogous nonresonant ion-atom systems. 



A. Implications to experimental data 

In the experiments by Kohl and coworkers 0] the ki- 
netic energy of a single Yb + ion immersed in an ultracold 
Rb enesemble was varied by adding excess micromotion 
energy after displacement of an ion from the center of a 
trap. The binary-collision ion-loss rate coefficient deter- 
mined in this way does not correspond to a conventional 
thermal rate constant for the Maxwell collision energy 
distribution at certain temperature. For this reason, we 




co (cm ) 



FIG. 6: (color online) Frequency-resolved radiative cross sec- 
tions normalized to the total cross sections, FGR calculations 
for the SO model. Lines correspond to the collision energy 
of 0.01 cm - , shaded areas - 1 cm -1 . Vertical dotted lines 
mark the borders of free-bound spectrum set by the ground 
rovibrational level and dissociation limit of the X state. For 
RA, only the transitions with intensities exceeding 1% of the 
maximum one are shown. 

also used the effective energy-dependent rate coefficient 

R(E) = ^2E^a R AX (E) 

omitting negligible nonradiative CT contribution. This 
quantity derived from the OP calculations with the SO 
model is displayed in Fig. 7 within the energy range 
probed experimentally [1]. In agreement with the mea- 
surements and the trait of Langevin regime, the back- 
ground rate coefficient does not depend on a collision 
energy. To define the average value, we took an arith- 
metic mean over an uniform grid of 10000 points covering 
the relevant energy range from 0.15 to 3.25 cm" 1 . For 
the 174 Yb ion, it amounts to 2.9 x 10~ 14 cm 3 /s in good 
agreement with the measured mean of (4.0 ±0.3) x 10~ 14 
cm 3 /s [1]. The background is ca. 2.6 x 10 -14 cm 3 /s, so 
the resonances contribute around 10% to the total rate. 
The SR model gives the mean rate of R — 2.4 x 10~ 14 
cm 3 /s. 

Experiments with the 172 Yb + ions revealed a quite sig- 
nificant isotope effect, -R174/-R172 ~ 1-4, see Fig. 7. Using 
the OP approximation Q and assuming that the terms in 
the sum do not depend on the reduced mass, one arrives 
to small reverse scaling, R174/R172 cx (/im/Vm) 3 / 2 ~ 
0.994. It is supported by the fact that calculated back- 
ground rate coefficients for two isotopes are hardly distin- 
guishable. Nevertheless, the mean rate coefficients shown 
in Fig. 7 reproduce the reverse isotope effect. It originates 
from the resonant contribution and reflects slightly lower 
resonance density for the lighter isotope (careful inspec- 
tion of Fig. 5 indeed reveals the less congested structure) . 
Still, the theoretical rate coefficient ratio 1.03 is much 
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FIG. 7: (color online) The energy-dependent CT rate coef- 
ficient Hi 74(E) (the SO model, OP calculations). Solid and 
dashed (red) lines represent the mean values for collisions with 
the 174 Yb + and 172 Yb + ions, respectively. Shaded areas rep- 
resent the experimental error bars [3]. Rate coefficients are 
shown in the collision energy range studied experimentally. 



smaller than the measured. More recent experiments for 
the same system have indicated dramatic effect of the 
hyperfine structure [5] 169) . which likely contributes to 
isotope effect as well. 

Despite very reasonable agreement for the rate coeffi- 
cient, our results do not reconcile to the measured prod- 
uct distributions. Experimentally, the probability of Rb + 
ion production was found to be 35%, while the rest 65% 
correspond to the loss of charged particle from the trap 
HIE]. Our FGR calculations also give 30% for CE prob- 
ability equivalent to Rb + ion formation, but assign the 
remaining 70% to RA, formation of the stable (YbRb) + 
molecular ion not observed experimentally. Our calcu- 
lations do not identify the trap loss channel, because, 
according to Fig. 6, all the CE products have the kinetic 
energy less than 500 cm -1 at the trap depth of 1200 
cm -1 @]. Several works [SHU [70] emphasized efficient 
secondary collision processes for Ca + , Ba + and Rb + ions 
in contact with an ultracold Rb ensemble. This reason, 
however, was ruled out by Ratschbacher et al. [8], who 
failed to observe (YbRb) + ion even at very short con- 
tact times. We assume that a more plausible explanation 
can be efficient photodissociation of the molecular ion to 
highly energetic products by cooling or trapping optical 
fields. More quantitative estimation may be useful for 
bridging the gap between theory and experiment. 



B. Other nonresonant cold ion-atom collisions 

In many respects, neutral and singly ionized Yb is 
the 4/-analog of an alkaline-earth metal (M). Indeed, a 
lot of similarity can be traced out in the interactions 
of Rb atom with Yb+ and M+ ions. Due to mismatch 



in ionization potentials (IPs), the entrance M + ( 2 S) + 
Rb( 2 S) channel always lies above the ground M(*S) + 
Rb + ( 1 S) one, so the location and symmetry of excited 
M* + Rb+( 1 S) CT channels determine the main qualita- 
tive difference. Binding of ns 2 electrons decreases with 
n from Be to Ba as manifested in IP, ns 2 — > nsnp and 
ns 2 — > (n — l)dns promotion energies [46) . By these pa- 
rameters, Yb most closely resembles Ca. However, in 
contrast to Yb, the excited Ca*( 3 P°) + Rb CT chan- 
nel is open and lying slightly below the entrance [7j [30] . 
Moszynski and coworkers [3TJ studied theoretically the 
Ba + + Rb system, in which excited CT channels are 
energetically closed, see also Ref. [71] . However, these 
states correlate to Ba*( 3 D) + Rb + ( 1 S) asymptote and 
exhibit more complicated interaction with the entrance 
y4 1 S+ and a 3 E+ states. In particular, the PEC of the A 
state has a double-well shape even in the nonrelativistic 
case. It should be expected that the Sr + + Rb system 
reveals more similarity to the Yb + + Rb one, but we are 
not aware of any work on it. 

It appears that Yb + + Rb represents one of the sim- 
plest systems from the viewpoint of excited state effects. 
The use of lighter alkali partners may further simplify the 
dynamics because increasing IP of the neutral will push 
the entrance channel down in energy. In the extreme case 
of Li, it lies only 7000 cm -1 above the ground and almost 
10000 cm -1 below the excited CT asymptote. Collisions 
with Li should be of special interest for cooling a trapped 
ion below the Langevin regime, as discussed by Cetina et 
al. [25]. 

Rellergert et al. [6] investigated cold Yb + + Ca col- 
lisions, in which only two states of the 2 E + symme- 
try separated asymptotically by 1140 cm -1 are involved. 
Avoided crossing at long range presents an interesting 
feature of this system that facilitate nonadiabatic CT 
pathway, see below. Similar picture should be expected 
for Yb+ + Sr, whereas for Yb+ + Ba the CT to ex- 
cited Ba + ion becomes possible. Among the nonresonant 
alkaline-earth systems, Ba + + Ca has been studied by 
Sullivan et al. [26J. Here the CT process is endothermic 
and occurs from excited states of the ion. 

It is also instructive to bring together scarce, mostly 
theoretical, data on the efficiency of distinct CT path- 
ways. To do so on the same footing, we presented cl 
coefficients - the ratios of the rate coefficient (cross sec- 
tion) to its Langevin inelastic value. 

Our estimation for nonradiative CT rate is vanishingly 
small, £l « 10~ 17 . Much larger efficiencies, £l ~ 5xl0 -6 
and 10 -3 , were reported for Yb + + Ca [5] and Ca + + 
Rb [5H1 133], respectively. In the former case, nonadi- 
abatic transitions should be dramatically enhanced by 
above mentioned avoided crossing. The latter coefficient 
in fact refers to the transition to the open excited CT 
channel, while the efficiency of the direct A — X CT is 
of order of 10~ 6 [30]. The measured e L « 10~ 3 for Ca+ 
+ Rb [10] confirms the nonradiative CT to excited chan- 
nels. Indeed, the efficiency of the direct radiative CT was 
found to be by an order of magnitude lower, ej, k 7x 10 -5 
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ma- 

For radiative CT we got €l ~ 1.5 x 10 5 , in good 
correspondence with Ca + + Rb and Ba + + Rb collisions 
(theory e L m 4 x 10" 5 [3T] and 5 x 10" 6 [71]; experiment 
10~ 4 - 10" 3 [5] and < 2 x 10" 4 [71]). Calculations for 
Yb+ + Ca provided ej, s» 0.03, also in agreement with 
the measured rate [5] . The lowest radiative CT efficiency 
£l ~ 10~ 7 was calculated for the Ca + + Na collisions 

[nun]. 

Even less is known on the product distributions. The 
molecular ion was observed in the Ca + + Rb system to- 
gether with RbJ ion formed in the secondary collisions 
□ . Studies of the Ba+ + Rb [7T| and Yb+ + Ca [5] 
collisions revealed the same discrepancy between experi- 
ment and theory as we met here. In the former case, the 
(BaRb) + ion constitutes 30% of the ionic product, while 
the theory designates it as the barely dominant product 
of the radiative CT. An upper bound for (YbCa) + forma- 
tion probability was estimated experimentally as 0.02%, 
whereas according to FGR calculations radiative associ- 
ation to molecular ion contributes 50%. It was specu- 
lated that the secondary RA processes may deplete the 
diatomic ion to heavier species like (Ca2Yb) + [Jjj. In our 
case the measurements failed to detect charged particles 
different from atomic ions and we proposed alternative 
explanation - photodissociation of a molecular ion by 
cooling and trapping optical field. Finally, calculations 
on Ca + + Na collisions showed even stronger predomi- 
nance of the radiative molecular ion formation: RA to 
CE branching ratio was found to be around 20:1 [T3lll7j . 

The above overview clearly reveals a variety of the 
CT mechanisms in cold ion-atom collisions and orders- 
of-magnitude variations of their rates. The obvious rea- 
son is the strong dependence on tiny electronic proper- 
ties of a system: a relative potential energy gap between 
initial and final CT channels, transition dipole moments, 
nonadiabatic couplings. It can be inferred that the radia- 
tive pathway generally dominates the direct CT process 
and leads mainly to the formation of molecular ion. The 
nonradiative pathway can be competitive if CT occurs to 
excited states of the products lying close from below to 
the entrance. The CT processes involving excited ions or 
atoms [BHSJ [HI HJO] may well exhibit a complex interplay 
of the distinct pathways. 

VI. CONCLUSIONS 

High-level scalar-relativistic ab initio calculations on 
the lowest excited electronic states of the (YbRb) + ion 
identified the direct coupling between the entrance A 1 Y J + 
and exit X 1 Tj + states as the only pathway for charge 
transfer in the cold Yb + ( 2 S) + Rb( 2 S) collisions. Spin- 
orbit coupling leaves the charge transfer through the en- 



trance al( 3 E + ) state inefficient, but modifies the poten- 
tial energy curve of the A state to the double- well shape 
due to avoided crossing with the excited state correlating 
to the closed Yb*( 3 P°) + Rb + charge-transfer asymp- 
tote. 

Quantum scattering calculations showed that the non- 
radiative charge transfer induced by radial and angular 
nonadiabatic coupling matrix elements has a negligible 
transition probability. Charge transfer proceeds radia- 
tively with the effective rate coefficient ca. 3 x 10~ 14 
cm 3 /s, that is, five orders of magnitude lower than the 
Langevin capture rate. The calculated radiative cross 
sections bear all traits of the Langevin regime includ- 
ing the rich structure associated with centrifugal barrier 
tunneling (orbiting) resonances. The short-range well of 
the entrance A0 + state has a weak effect on the charge- 
transfer efficiency due to decline of the transition mo- 
ment, but further complicates the resonance structure. 

Calculated radiative charge-transfer rate coefficients 
agree well with the values of (4.0 ± 0.3) x 10 -14 and 
(2.8 ± 0.3) x 10~ 14 cm 3 /s measured for 174 Yb+ and 
172 j onS; respectively @]. The observed isotope ef- 
fect is opposite to that expected from the mass factor for 
a radiative charge-transfer rate. In our calculations, it is 
reproduced correctly and originates from the resonance 
contribution, which increases with the mass of the sys- 
tem. The magnitude of the effect, however, was found to 
be too small to fully explain the observations. 

The calculations estimated the probability to find the 
product ion Rb + in a trap as ca. 30% in agreement with 
35% obtained experimentally [4, . The rest of the events 
correspond to the formation of the bound (YbRb) + 
molecular ion, in sharp contrast to observed loss of the 
charged particle. Efficient molecular ion photodissocia- 
tion to highly energetic products by cooling and trapping 
fields may explain this disagreement. 

Present results are in line with the previous studies of 
similar processes involving alkaline-earth ions and atoms. 
They suggest that the radiative pathway dominates the 
ground-state charge transfer having a mean efficiency 
around 10 -5 of the Langevin rate and the strong propen- 
sity to the molecular ion formation. 
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